Hysteresis and avalanches in the T = random-field Ising model with 2-spin-flip 

dynamics. 



Eduard ViveiQ 

Dept. d'Estructura i Constituents de la Materia, Universitat de Barcelona 
Diagonal 647, Facultat de Fisica, 08028 Barcelona, Catalonia 

Martin Luc Rosinberg and Gilles Tarjus 
Laboratoire de Physique Theorique des Liquides, Universite Pierre et Marie Curie, 4 Place Jussieu, 75252 Paris, France 

(Dated: February 2, 2008) 

We study the non-equilibrium behavior of the three-dimensional Gaussian random-field Ising 
model at T = in the presence of a uniform external field using a 2-spin-flip dynamics. The 
deterministic, history-dependent evolution of the system is compared with the one obtained with the 
standard 1-spin-flip dynamics used in previous studies of the model. The change in the dynamics 
yields a significant suppression of coercivity, but the distribution of avalanches (in number and 
size) stays remarkably similar, except for the largest ones that are responsible for the jump in the 
saturation magnetization curve at low disorder in the thermodynamic limit. By performing a finite- 
size scaling study, we find strong evidence that the change in the dynamics does not modify the 
universality class of the disorder-induced phase transition. 
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I. INTRODUCTION 

The non-equilibrium random field Ising model (RFIM) 
was introduced by Sethna et al.^ as a model for the 
Barkhausen effect in ferromagnets and more generally as 
a prototype for many experimental systems that show 
hysteretic and jerky behavior when driven by an exter- 
nal force. Because of the presence of disorder, these 
systems have a "complex free energy landscape" with a 
multitude of local minima (or metastable states) sepa- 
rated by sizeable barriers, which makes thermally acti- 
vated processes essentially irrelevant at low enough tem- 
perature (the lifetime of metastable states may then be 
considered as infinite). As a consequence, these systems 
remain far from equilibrium on the experimental time 
scales (even when the driving rate goes to zero) and 
their response to the external force is made of a series 
of jumps (avalanches) between neighboring metastable 
states. This type of behavior is very well modeled by the 
ferromagnetic RFIM with a zero-temperature single-spin- 
flip dynamics in which a spin flips only if this lowers its 
energy. The local character of the energy minimization is 
then at the origin of irreversibility. With this dynamics, 
the RFIM satisfies the property of return-point memory 
(or "wiping out" effect) which is a feature observed in 
several experimental systems with good approximation. 
Moreover, in dimension d > 3, the model is known to 
exhibit an out-of-equilibrium phase transition between a 
strong-disorder regime where the magnetization hystere- 
sis loop is smooth on the macroscopic scale and a weak- 
disorder one where it has a discontinuous jump. Such a 
transition has been observed in thin Co/CoO filmsH and 
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Cu-Al-Mn alloysQ, and it has been recently suggested^ 
that it may also be associated to the change in the ad- 
sorption behavior of '*Hc in dilute silica aerog els0. The 
two regimes, strong and weak disorder, are separated by 
a critical point characterized by universal exponents and 
scaling laws which have been extensively studied by ana- 
lytical and numerical methods|^0]. In particular, much 
effort has been recently devoted to analyze the number 
of avalanches, their size, and their geometrical properties 
above, below, and at criticality|^ These results pro- 
vide a comprehensive, though rather complex scenario 
for the phase diagram of the non-equilibrium RFIM with 
a metastable dynamics in the thermodynamic limit. A 
recent discussion of the relevance of the model to the de- 
scription of the Barkhausen effect in real magnets can be 
found in Refs.[l3,[l3. 

An issue that so far has not been studied is the ro- 
bustness of this theoretical description with respect to a 
change in the dynamics (in the literature on the RFIM, 
it is implicitely taken for granted that this should not 
matter). The single-spin-flip dynamics, however, is not 
the unique (and may be not the best) way of simulat- 
ing hysteretic dynamical processes in actual systems. It 
is clear for instance that the hysteresis loop will shrink 
if the dynamics allows for a better equilibration of the 
system by employing multiple-spin flips. Then, what 
will be the avalanche properties ? Will there still be 
a phase transition ? If so, will the critical behavior be 
the same as with the single- spin- flip dynamics ? There is 
in fact the intriguing possibility, supported by numerical 
simulations and analytical arguments 1^ Il2l| . that the 
non-equilibrium and equilibrium transitions of the T ~ 
RFIM belong to the same universality class, even if crit- 
icality occurs in zero external field at equilibrium and at 
a non-zero coercive field in the irreversible evolution [T^. 
Since the ground state is stable with respect to the flip of 
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an arbitrary (finite) number of spins, this may indicate 
that the disorder-induced transition has a universal char- 
acter at criticahty which does not depend on the specific 
choice of the dynamics 0. 

In order to shed some light on this issue and check 
the robustness of the transition, we study here the non- 
equilibrium r = RFIM with a 2-spin-flip dynamics. We 
compare the results with those obtained with the stan- 
dard 1-spin-flip dynamics, in particular those concerning 
the number and size of the avalanches. We first show in 
section II that one can indeed define a 2-spin-flip algo- 
rithm that yields a deterministic evolution of the system 
with the external field. In particular, the dynamics sat- 
isfies the "abelian" property which guarantees that the 
same 2-spin-flip stable configuration is attained, what- 
ever the order in which the unstable spins are relaxed 
during an avalanche. It also has the property of return- 
point memory. In section III we present the results of 
our numerical simulations on 3-d lattices and compare 
them to the behavior of the same system with a 1-spin- 
flip dynamics. The hysteresis loops arc signiflcantly re- 
duced when allowing 2-spin flips, but the main features, 
including the presence of a disorder-induced transition 
with an associated critical point, are not significantly al- 
tered. We then perform in section IV a finite-size scaling 
analysis, which allows for a determination of the critical 
properties. Wc find that, within statistical uncertainty, 
the exponents and scaling functions are identical to those 
obtained with the standard 1-spin-flip dynamics. The 
main conclusions of the study are reported in section V. 



II. MODEL AND DYNAMICS 

The model is defined on a cubic lattice of linear size 
L with periodic boundary conditions. On each site (z = 
1, . . . , = L^) there is an Ising spin variable (5*^ = ±1). 
The Hamiltonian is 

<i,j> i i 

where the first sum extends over all distinct pairs of 
nearest- neighbor (n.n), H is the external applied field, 
and hi are quenched random fields drawn independently 
from a Gaussian distribution with zero mean and stan- 
dard deviation a. We arc interested in studying the se- 
quence of states along irreversible paths at T = when 
the system is driven by the external field H (in the adi- 
abatic limit corresponding to a vanishingly small rate 
of change of the external field). For this purpose, the 
Hamiltonian must be supplemented by some dynamical 
rules. 

The standard 1-spin-flip dynamics used in previous 
studies consists in minimizing the energy of each spin, 

Ti-i — —Sifi (2) 
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FIG. 1: Stability diagrams showing the state with minimum 
local energy according to the values of the fields created by 
the neighborhood (defined in the text), (a) corresponds to a 
single spin i and (b) to a pair of neighboring spins i,j. 

where 

f, = jY,Sj+h, + H (3) 

is the net field at site i (the summation in Eq. (3) is over 
the z n.n. of i). From the above expression, it is clear 
that the minimization of Tii is obtained by aligning each 
spin with its local field. Si = sign(/i), as represented 
schematically in Fig. Q^a) . This provides a stability cri- 
terium for any state with respect ot this 1-spin-flip dy- 
namics. This dynamical rule may be implemented b y an 
algorithm that propagates one avalanche at a time[Ij]. 
Starting from a stable configuration, one increases (or de- 
creases) the external field until the local field fi at some 
site i becomes zero (this corresponds to the vanishing 
of the local minimum in which the system was trapped). 
The spin Si (which is uniquely defined because the distri- 
bution of the random fields is continuous) is then flipped, 
which in turn may cause neighboring spins to become 
unstable and thus initiate an avalanche. The avalanche 
stops when a new mctastable state is reached. The ex- 
ternal field is then changed again, and so on. When the 
spins that become unstable during the avalanche are se- 
quentially reversed (e.g., by increasing i from 1 to N), it 
is of course crucial that the final state does not depend 
on the sequential order. Thanks to the ferromagnetic na- 
ture of the couplings, this is indeed the case as a result 
of the so-called " no-p assing" and "abelian" properties of 
this dynamics^, Il5| . Moreover, the same state is also 
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reached when all unstable spins arc flipped in parallel, 
which allows to measure the "time" it takes an avalanche 
to occur 
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FIG. 2: Evolution of the enthalpy and the magnetization (per 
spin) for a system of size L = 30 and a — 2.5 as the external 
field is changed from H — +oo to H — —2. The trajectories 
corresponding to sequential and parallel updating algorithms 
are also compared with the adiabatic one. Points indicate the 
intermediate unstable states. The inset shows a closeup of 
the region around the endpoint. Notice that the end-point is 
the same in all cases. 

Setting the rules for a 2-spin-flip dynamics is rather 
straightforward. By definition, 2-spin-flip stable states 
are spin configurations whose energy (defined solely by 
the Hamiltonian) cannot be lowered by the flip of one 
or two spins (clearly, new features are only introduced 
when these two spins are n.n.). The (local) energy to be 
minimized is thus the one associated with a pair ij of 
n.n. spins, 



fi fj J 



where 



n = j 



H 



(4) 



(5) 



is the fleld experienced by Si without the influence of the 
neighbor Sj (the summation in Eq. (5) is over the n.n. 
of i excepting j). One can then think of the dynamics 
as made of single-spin flips and "irrcducibly cooperative" 
2-spin flips. As pictured in Fig. H^b), a single-spin flip 
occurs whenever the net field on Si, fi = f- + JSj, or 
on Sj, fj = fj + J Si, changes sign. This corresponds 
to the changes ii^Ti, ii^it, TT^Ti, TT^iT in the 
diagram. An irreducibly cooperative 2-spin fiip involves 
a n.n. pair of spins with the same sign (ff or jj,) that 
cannot fiip individually (i.e. without the simultaneous 
flip of the neighbor). This occurs whenever the net field 



on the pair of aligned spins. = -|- /j, changes sign 
in the region of the diagram where —J<fl<J and 
— J < fj < J. These two last conditions come from the 
fact that the spins cannot individually fiip (hence neither 
/; nor fj changes sign before the cooperative flip of the 
pair) and the condition that once a pair has flipped, none 
of the spins can flip back individually. 

The corresponding algorithm is a simple extension of 
the one described above for the 1-spin-flip dynamics. 
Starting from a 2-spin-flip stable configuration, the ex- 
ternal field is varied until one finds a pair of spins that 
becomes marginally stable: the representative point of 
this pair in the diagram of Fig. ^h) leaves the region 
where it was originally (associated with ft, t|, |t, or H) 
and. depending on the border which is first attained, only 
one spin flips or the two spins flip simultaneously [ig. 

It is not hard to show that this new dynamics obeys the 
same properties as the 1-spin-flip dynamics, in particular 
the crucial abelian property. This is again a consequence 
of the ferromagnetic nature of the interactions. One only 
needs to note that the state of the system can be repre- 
sented by a set of zN/2 points (corresponding to all the 
distinct n.n. pairs) in the diagram of Fig. ^b). As the ex- 
ternal field is monotonously increased (resp. decreased), 
the local fields can only increase (resp. decrease), the 
spins can only fiip up (resp. down), and the points can 
only move up and right (resp. down and left) in the di- 
agram. By slightly modifying the arguments of Ref . Q , 
one then can prove the no-passing rule, the abelian prop- 
erty and the existence of return-point memory. 
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FIG. 3: Major hysteresis loop and an internal loop (—0.85 < 
H < 0.8) obtained with the 2-spin-flip dynamics for a system 
of size L = 30 and a = 2.5. The inset shows the details around 
the point at H = —0.85 revealing the property of return-point 
memory. 

Instead of paraphrasing the demonstrations given in 
Ref. , we choose here to illustrate these properties by a 
numerical example. (We have also performed numerical 
tests in many situations and found no violations of these 
properties.) The evolution of a system with size L = 30 
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and a = 2.5 is shown in Fig.|21where the energy per spin, 
e = Ti./N, is plotted as a function of the magnetization 
m = (strictly speaking, e is the enthalpy). The 

external field H is varied from a very large initial value 
where all spins are up to the final value H = —2 (here 
and after, J is taken as the energy unit). Two of the 
curves display the sequence of unstable states that arc 
obtained after a sudden change of the external field using 
either a sequential or a parallel updating algorithm. We 
also show the metastable evolution corresponding to the 
adiabatic driving (with sequential updating) along the 
hysteresis loop. In all cases the final state is the same. 
This is true even when the intermediate states are dis- 
tinct, for instance when the spins are chosen sequentially 
in a different order. 

The property of return-point memory property is il- 
lustrated in Fig. 121 again for a system with size i = 30 
and a = 2.5. The minor loop is obtained by reversing 
the evolution of the external field first in the decreasing 
branch at H = —0.85 and then a± H ~ -f-O.SO. As can 
be seen from the inset, the internal loop closes before the 
return point, so that the evolution follows that of the 
major loop for a small region oi H > —0.85. 



III. NUMERICAL SIMULATIONS 
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FIG. 4: Magnetization curves obtained with the 1-spin and 
2-spin-flip dynamics in a sample of size L = 30 for (a) cr = 2 
and (b) (7 = 3: dotted-dashed lines correspond to the 1-spin- 
fiip dynamics and dashed lines to the 2-spin-flip dynamics. In 
addition, the ground-state magnetization is shown as contin- 
uous lines. 



We now present the results of numerical simulations 
performed on 3-dimcnsional cubic lattices using either 



the 1-spin-flip or 2-spin-flip dynamics with sequential up- 
dating. Averaged quantities were obtained with statistics 
over 10'^ to 10^ different realizations of the random field 
distribution and system sizes ranging from L = 8 up to 
L = 48. As emphasized in Ref.ja], in order to describe 
properly avalanche properties (especially the "spanning" 
avalanches), it is more important to perform averages 
over many disorder realizations than to simulate very 
large system sizes. 



A. Hysteresis loops 
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FIG. 5: Enthalpy per spin as a function of H corresponding 
to the loops in Fig. 5 (for clarity, only the ascending branches 
are shown). The continuous lines represent the ground-state 
behavior. 

Fig. 01 shows the hysteresis loops obtained in a single 
sample for two different values of the disorder a. For 
comparison, wc also display the magnetization curves 
obtained with the 1-spin-fiip dynamics and with the al- 
gorithm of Ref.^3| which gives the exact ground-state 
(equilibrium) magnetization. The corresponding behav- 
ior of the enthalpy per spin along the ascending branches 
of the loops is reported in Fig. jsl (note in passing that 
the ground-state enthalpy does not show any disconti- 
nuity as the external field is variedflTj'l. As could be 
expected, the main effect of the new dynamics is to re- 
duce the size of the hysteresis loops. Specifically, the 
cocrcivity (i.e., the magnitude of the external field for 
which the magnetization is equal to zero) is decreased by 
more than 30% when allowing pairs of spins to flip to- 
gether. Accordingly, the enthalpy difference between the 
ground state and the metastable states which are visited 
along the loops is also reduced. Nevertheless, the loops 
display the same key feature, that is a change from a 
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discontinuous to a continuous behavior as the disorder 
is increased. This suggests that there is also an out-of- 
equilibrium disorder-induced phase transition under the 
2-spin-flip dynamics, with a critical value of a at which 
the discontinuity appears for the first time in the ther- 
modynamic limit. 

It is worth pointing out that the new dynamical rules 
allows the system to effectively overcome energy barriers 
of magnitude up to A_E = 2J. Indeed, the difference 
between the two dynamics shows up when a pair of n.n. 
spins with same sign can cooperatively flip, say from J, J, 
to tt when the external field is increased, whereas each of 
its spin cannot individually flip. This means that along 

1- spin-flip paths, the system has now been able to by- 
pass the higher-energy states, cither f J. or |f . By using 
Eq. (4), it is easy to show that the relevant barrier height 
associated with this process is at most 2 J. Since coopera- 
tive flips of more than 2 spins do not occur with the cho- 
sen dynamics and the system's trajectory otherwise go 
through states of decreasing energy, one concludes that 

= 2 J is the maximum barrier height that the system 
may overcome when passing from the 1-spin-flip to the 

2- spin-flip dynamics. 



B. Avalanches 

As shown in recent studies [^Q, a good characteriza- 
tion of the disorder-induced critical point can be reached 
by analyzing the number and size distribution of the 
magnetization jumps (avalanches) that compose the hys- 
teresis loops in finite systems. For that purpose, it is 
necessary to classify the avalanches in several categories, 
according to their behavior as the system size L is in- 
creased. One first has to distinguish whether or not an 
avalanche spans the system from one side to the other, 
in 1, 2 or 3 spatial directions (indicated in the following 
by the index a). For each individual avalanche, this is 
a property that can be easily detected during the sim- 
ulation. Avalanches are thus classified as being non- 
spanning [a =ns), Id-spanning (a = 1), 2d-spanning 
(a = 2), or 3d-spanning (a = 3). 

Fig. shows the number of Id, 2d and 3d-spanning 
avalanches recorded along the descending branch of the 
hysteresis loops as a function of a. The data, averaged 
over disorder, correspond to a system of size L = 24. It 
can be seen that the behavior of the three quantities is 
completely equivalent under the two dynamics. The only 
difference is a shift towards larger values of a when the 
2-spin-flip dynamics is used. The same shift is also found 
for all studied system sizes. This is a first indication that 
(j^c^ > ai^"^ , as will be confirmed by the finite-size scaling 
analysis presented in the next section. In the case of the 
1-spin-flip dynamics, a detailed analysis was performed in 
Refs.|3,l3, revealing the scenario that occurs in the ther- 
modynamic limit and that is already suggested by the 
data shown in Fig. El when L — > oo, Ni{(j) and N2{<j) 
are expected to display a (5-singularity at Uc and N^{a) 



a step-like behavior. It was shown, moreover, that there 
are two types of 3d-spanning avalanches, subcritical and 
critical, which scale with different exponents. The former 
are responsible for the discontinuity in the magnetiza- 
tion curve in the thermodynamic limit (there is only one, 
compact, subcritical avalanche for cr < CTc) whereas the 
latter only exist at CTc (hence the additional ^-singularity 
at the edge of the step function whose signature is al- 
ready visible in Fig. IBJ. In a flnite system, however, all 
kinds of avalanches may exist close enough to the critical 
point, and it is quite difficult to discriminate subcritical 
from critical avalanches. In Ref.|^, an elaborate analysis 
was needed to show that these avalanches have different 
fractal dimensions at criticality. This study is impossible 
here because of the complexity of the 2-spin-fiip algo- 
rithm that forbids the use of large systems with good 
enough statistics. Therefore, in the following, we shall 
not distinguish between these subcritical and critical 3d- 
spanning avalanches. 
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FIG. 6: Average numbers of Id, 2d, and 3d-spanning 
avalanches as a function of cr in a system of size L = 24. 
Continuous and dashed lines are guides for the eye and corre- 
spond to the 1-spin-flip and 2-spin-flip dynamics, respectively. 

Fig. 13 shows the (normalized) avalanche size distribu- 
tions Da{s;a,L) obtained along one branch of the hys- 
teresis loop for three different values of a (for clarity, 
the results obtained with the 1-spin-fiip dynamics are 
represented by continuous lines). Surprisingly, one can 
see that the distribution of non-spanning avalanches in 
Figs. CJa) to [Tfc) is almost unaltered by the change in 
the dynamics. The only small difference (barely visi- 
ble on the figure) induced by the 2-spin-fiip dynamics is 
that there are a little less avalanches of size s = 1 and 
a little more avalanches of size s = 2, but the rest of 
the distribution is almost the same. In particular, with 
both dynamics, the expected power-law behavior of the 
distribution will be characterized by the same exponent 
Teff w 2.0 in the thermodynamic limit|^. 

The size distributions Di and D2 of the Id and 2d- 
spanning avalanches shown in Figs. [TJd) to Efg) also 
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FIG. 7: Avalanche size distributions for a = 1.60 (a,h), 2.25 (b,d,f,i), and 2.80 (c,e,g,j). The open symbols and the continuous 
lines correspond to the 2-spin-flip and 1-spin-flip dynamics, respectively. The first row shows the distribution Dns of non- 
spanning avalanches and the other rows show the distributions D\ , D2 and Ds of the Id, 2d and 3d-spanning avalanches, 
respectively. All data have been obtained in a system of size L = 24. In (b), the dashed line indicates the expected power-law 
behavior of Dns at criticality with exponent Te// = 2.0. 



appear to be identical with the two dynamics, at least 
within statistical error bars. (Note that these avalanches 
do not exist for cr = 1.6 because this value is much lower 
than (Tc for both dynamics.) The only visible differ- 
ences between the two dynamics occur in D^, the size 
distribution of the 3d-spanning avalanches. Specifically, 
for (T = 1.60 and 2.25 (i.e. below CTc and very close 
to (Tc, respectively), the large 3d-spanning avalanches 
tend to be shifted to even larger sizes. According to 
Refs. 0, 0, these avalanches are probably subcritical 
spanning avalanches and their average size is thus a mea- 
sure of the order parameter. Therefore, this result is 
another indication that ai^"^ > a^^K In contrast, for 
a = 2.80 (which is clearly above dc), the distribution D3 
is not affected by the dynamics (Fig[7{j)): in this case, 
one expects to detect only critical 3d-spanning avalanches 
in a finite system. 



IV. FINITE-SIZE SCALING ANALYSIS 

(2) 

As already mentioned, a precise determination of CTc 
and of some of the critical exponents requires a detailed 
finite-size scaling analysis which, unfortunately, is not 
posssible with the present algorithm. What can be done, 
however, is to check whether the set of exponents found 
with the 1-spin-flip dynamics can also be used to scale 
the new 2-spin-flip data. We first consider iVi and N2, 
the numbers of Id and 2d-spanning avalanches, which are 
quantities that have a simple scaling behavior. Following 
previous work|^|^, we assume the forms 

iVi (cr, L) = L'^Ni (^2^^/") (6) 

7V2(a,L) = L«7V2 (u2^'/") (7) 

where 9 and v are critical exponents, iVi and N2 are scal- 
ing functions, and U2 is a scaling variable that measures 
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the distance to the critical point. It is defined as 



J2) 
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A 



(8) 



where the parameter A accounts for a second order cor- 
rection that plays a role when the studied systems are 
not very large, as is the case here. In Ref.Q, the best 
choice for the collapse of the scaling plots was obtained 
with V = 1.2, 9 = 0.1, and A = —0.2, values that we keep 

(2) 

here. The only free parameter is thus CTc ■ 

As shown in Fig. |S1 a very good collapse of the new 

(2) 

data can be obtained with CTc = 2.25. Taking into ac- 
count the quality of the plots, there is an uncertainty of 
±0.01 on this value, but, clearly, the value ctc^'' = 2.21 ob- 
tained with the 1-spin-flip dynamicsQ can be discarded. 
We have also checked that this conclusion is not modified 
when the non-universal parameter A is allowed to vary. 
These scaling plots thus indicate that the 2-spin-flip dy- 
namics only induces a shift of the critical value of the 
disorder but does not change the universality class of the 
transition. 

A second check of this result may be obtained by mea- 
suring the average field (i?3(cr)) at which the 3d-spanning 
avalanches occur. This quantity was studied with the 
1-spin-flip dynamics in Ref. 0, which allowed to map 
out the first-order line (corresponding to the macroscopic 
jump in the magnetization) in the diagram H ~ a. The 
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FIG. 9: (a) Average field {Hz{(y, L)) at which the Sd-spanning 
avalanches occur as a function of a for different system sizes, 
(b) Scaling plot of the data below (j'^^ according to Eq. 1^1 
with JfP' = -0.885, B' = 0.25 and = 1.5. 



dependence of {H:i{a)) with the system size is shown in 
Fig. EJa). Note that the transition line extends above 
in a finite system. However, the end-point, beyond which 
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no 3d-spanning avalanches are found (with the present 
samphng of disorder realizations), becomes closer and 
closer to the critical point {He, cTc) as L is increased. 
According to Ref. j^l , this set of curves should scale as 



J2)' 



.(2) 



(9) 

(2) 

where He is the critical field, B' a non-universal tilting 
constant, a critical exponent, and the corresponding 
scaling function. Strictly speaking, the scaling should be 
done separately for the average fields (H3-) and (i?3c) at 
which the subcritical and critical 3d-spanning avalanches 
occur0. Indeed, the number of these avalanches scales 
differently with L. However we expect the lack of scaling 
to have a (small) effect only in the region U2L^^^ — 1 
where the two kinds of avalanches coexist 0- 

In Ref. , the best collapse of the 1-spin-flip data was 

obtained with ^ = 1.5, B' = 0.25, and H^^^ = -1.425. 
Here, we keep the same values for /i and B' , set dc = 
2.25, and consider the critical field as the only free pa- 
rameter. As shown in Fig. Elb), a very good collapse 

(2'] 

can be obtained with He — —0.885 (for clarity, we only 

(2) 

present the collapse for a < <Tc )■ Again, this result is 
consistent with the assumption that the new dynamics 
does not change the universality class of the transition. 
On the other hand, the significant decrease in the criti- 
cal field (in absolute value) is in line with the decrease in 
coercivity illustrated by Fig. 4. 

Finally, it is interesting to study the influence of the 
new dynamics on the finite-size scaling functions. In 
Fig. ^1 we compare the scaling collapses of the num- 
ber of Id-spanning avalanches, A^i (cr,i), obtained with 
the two dynamics (for the 2-spin-flip dynamics, this is 
the same curve as in Fig. 8). As can be seen, the agree- 
ment between the two curves is quite remarkable. Even 
the deviations (around U2L^^'^ ~ 0) from the Gaussian 
fit proposed in Ref.|^ are the same. A similar analysis 
for N2{(J,L) shows the same agreement. This is again a 
strong indication for universalit y, th at goes beyond the 
equality of the critical exponents [iq. 



V. CONCLUSION 

We have shown in this paper that the non-equilibrium 
behavior of the 3-d RFIM at zero temperature is not 
qualitatively altered when going from the standard 1- 
spin-flip metastable dynamics to a 2-spin-flip metastable 
dynamics. The coercivity (that is essentially the width 
of the hysteresis loop associated with the evolution of 
the magnetization with the driving field) is significantly 
reduced by allowing 2-spin fiips, but the main features 
of the hysteretic behavior, including the presence of a 
disorder-induced non-equilibrium transition and the dis- 
tribution of avalanches, remain similar. By using a finite- 
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FIG. 10: Comparison of the scaling plots of the number of 
Id-spanning avalanches obtained with the 1-spin-flip (sym- 
bols with error bars) and 2-spin-flip (symbols without error 
bars) dynamics. The two collapses only differ by the value 
of the critical disorder (cr?^ = 2.21 and af^ = 2.25). The 
continuous line is the Gaussian fit proposed in Ref. 



size scaling analysis, focused on the number and size of 
the avalanches, we have furthermore provided strong ev- 
idence that the critical behavior (exponents and scaling 
functions) obtained with the 2-spin-flip dynamics is in 
the same universality class as that obtained with the 1- 
spin-flip dynamics. 

Changing the dynamical rules used to study the evo- 
lution of the model, as done here, helps adressing several 
important questions. The first one is the robustness of 
the hysteretic scenario provided by the T = RFIM with 
the standard 1-spin-flip dynamics. As discussed in intro- 
duction, a basic assumption underlying the theoretical 
description is that the system gets trapped in metastable 
states on the experimental time scale and can only es- 
cape when a change in the external field makes the rele- 
vant state looses its stability. However, real systems are 
not at r = and some partial, local equilibration, due 
for instance to thermally activated processes, may take 
place even though the system remains far from equilib- 
rium on the experimental time scale. Introducing coop- 
erative 2-spin flips is a simple way to check the effect of 
partial equilibration processes. Our results clearly point 
towards the robustness of the whole theoretical picture 
drawn from previous studies of the model using the T = 
1-spin-flip dynamics 0, H, H H . 

A second question concerns the relation between the 
disorder-induced critical properties observed in the non- 
equilibrium behavior of the RFIM at T = and the 
equilibrium critical behavior associated with the param- 
agnetic to ferromagnetic transition. There is numerical 
evidence that, at least to a good approximation, criti- 
cal exponents and scaling functions associated with the 
two kinds of criticality arc the same 
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tional, but less conclusive evidence is provided by exact, 
but mean-field-like results on the Bethe lattice [T^l and 
by perturbation theory near the upper critical dimension 
d = 60.) Our present findings suggest that a whole se- 
ries of T = metastable dynamics involving fc-spin flips 
lead to the same non-equilibrium criticality, with the crit- 
ical disorder strength increasing with k and the critical 
coercive field decreasing with k. Equilibrium behavior 
at T = as a function of the external field involves the 
system's ground state, i.e. a state stable to flips of any 
arbitrary finite number of spins. Based on the numerical 
closeness of the critical exponents and scaling functions, 
on the fact that the critical disorder strength satisfies 
a^i > CT?' > ai^^ (whereas H^" = < H?' < i/i^^), 
and on the similarity of the underlying physics at T = 0, 
it is tempting to speculate that the equilibrium behav- 
ior can be obtained as the limit of a series of fc-spin-flip 



metastable dynamics with increasing k, with the critical 
properties of the whole series belonging to the same uni- 
versality class and governed by the same fixed point 0. 
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